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q \ Abstract 



We study the stochastic resonance (SR) of a flexible polymer chain 
crossing over a piecewise linear bistable potential. The dependence of 



Q\ \ signal to noise ratio SNR on noise intensity D, coupling constant k 

^^ ■ and polymer length N is studied via two state approximation. We 

find that the response of signal to the background noise strength is 



significant at optimum values of D opt , k op t and N opt which suggests 
<^ ' novel means of manipulating proteins or vesicles. Furthermore, the 

thermally activated barrier crossing rate r^ for the flexible polymer 
chain is studied. We find that the crossing rate r\. exhibits an optimal 
value at an optimal coupling constant k opt ; k opt decreases with N. As 
S^ ■ the chain length N increases, the escape rate for the center of mass r^ 

monotonously decreases. On the other hand, the crossing rate for the 
portion of polymer segment r s increases and saturates to a constant 
rate as ./V steps up. 
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1 Introduction 

Understanding the physics of thermally activated barrier crossing for systems 
consisting of one or more than one degree of freedoms [H[2],[3i[Ii[5j|6j[7J[8] is 
ubiquitous in many scientific disciplines. It has diverse physical applications 
and serves as a tool to understand stochastic paradigms such as stochastic 
resonance and thermal ratchets. Furthermore, for bistable systems, apply- 
ing fields such as an external load or nonuniform temperature breaks their 
symmetry along the reaction coordinate which may result in a unidirectional 
motion of particles. In the last few decades several model systems, which are 
far from equilibrium, have been proposed [HI HO] • Especially when the exter- 
nal load (force) varies in time, coordination of the noise with time varying 
force may lead to the phenomenon of stochastic resonance (SR) [TTl IT2] as 
long as the system is exposed to weak sinusoidal signals. 

The phenomenon of stochastic resonance (SR) is of interest these days 
because of its significant practical applications in a wide range of fields. SR 
depicts the phenomenon that systems, by utilizing the thermal background 
noise, enhance their performance when they are subjected to a small periodic 
signal. Since the pioneering work of Benzi. et.al [TT| . the idea of stochas- 
tic resonance has been broadened and implemented in many model systems 
[121 H3HHI HSlUni HH HSl HH I2D1 EI]- Recently the appearance of entropic 
stochastic resonance for a Brownian particle in a confined system has been 
reported in the work [22]. Unlike conventional energetic potential, the effec- 
tive potential of the confined systems may have an entropic nature due to the 
constrained regions. For tiny-scaled biological systems that frequently exist 
in a highly confined geometry, the entropic contribution to the effective po- 
tential is unavoidable and plays an important role in noise-induced resonant 
effects. 

Often biological and soft matter systems such as polymers and mem- 
branes are complex and consist of several components. Their flexibility and 
connectivity lead to a fascinating dynamics under certain time varying exter- 
nal fields and noisy environment. The resonance behavior of these systems 
relies not only on the strength of the background temperature, but also on 
their size, flexibility and shape of the potential. Thus, the interplay between 
the shape of the potential, flexibility and size of polymer (protein) plays 
a crucial role in the enhancement of signal to noise ratio SNR or spectral 
amplification rj as reported in the works [H)J dSJ H3, [23] • 

Earlier, Lindner et. al. considered linearly coupled damped bistable os- 



dilators [F5\. The dependence of SNR on the coupling strength k, number of 
oscillators N and noise intensity D was envisioned numerically. It has been 
shown that the SNR of the oscillators depicts a global maximum (SNR max ) 
at certain D and k op t for a given N. Latter, utilizing $4 field theory, March- 
esoni et.al [Hj independently checked the numerical results of Lindner et. al 
[TB] in large N limit. Recently, Dikshtein et.al [T7| considered a polymer in 
a symmetric bistable potential where the two end points of the polymer are 
restricted by motionless pining points in the perpendicular direction to the 
symmetric potential. It has been shown that the SNR is enhanced for wider 
(less deep potential). More recently we studied the stochastic resonance (SR) 
for a flexible polymer surmounting a bistable potential. Our analysis indi- 
cated that, due to the flexibility that can enhance crossing rate and change 
chain conformations at the barrier, the power amplification exhibits an opti- 
mal value at optimal chain lengths N opt and elastic constants k op t as well as 
at optimal noise strengths D |23j . 

The study of thermally activated escape rate of polymer has been also the 
subject of many studies (see for example [21] )• Since polymer is a flexible and 
an extended object with a finite length, its rate or SNR relies on its coupling 
constant, chain length, shape of the potential and initial conformation along 
the reaction coordinate in a complicated manner. Thus, more studies are 
needed to understand its complex dynamics. Most of the previous studies 
considered the center of mass motion. However, in this paper not only we 
examine the crossing rate for the center of mass motion r^ but also the rate 
for a portion of the polymer segment r s . We find that r^ < r s . 

The aim of this paper is to explore the crossing rate and stochastic res- 
onance of a flexible polymer chain in a piecewise linear bistable potential 
as a function of k, N and D by considering initially coiled chain conforma- 
tion. First we explore the escape rate of the chain as a function of the model 
parameters. We show that r^ monotonously decreases with N; the rate r s 
increases and saturates to a constant value when N steps up. Since the co- 
operation between the monomers increases with k, the rate increases as k 
goes up. At certain k opt , rj, attains an optimal value and further increases in 
k results in a lower rate as rigid polymer crosses the barrier at the expense 
of higher thermal kicks. 

In this paper, utilizing two state approximation, the dependence of SNR 
on D, k and N is examined. For globular polymer chain we show the response 
of signal to the background noise strength is significant at optimum values of 
D opt and N opt which suggests novel means of manipulating (such as efficient 



separation methods) not only for biopolymers, but also for proteins (vesicles) 
of different size. For the chain with a finite coupling constant k, the SNR 
exhibits an optimal value at an optimal k op t- The signal to noise ratio for 
the center of mass motion SNR^ exhibits an optimal value at an optimum 
N opt . On the other hand, the SNR S for the polymer segment monotonously 
increases with N. In addition, considering temperature dependent coupling 
constant (entropic chain) where the elastic constant is given by k = 3D /I 2 , 
we further confirm that SNR for entropic chain shows a broader and a higher 
peak than a chain with temperature independent coupling constant. Here, I 
designates the Kuhn segment length. 

At this point we stress that even though the model and its numerical 
approach is completely different from our previous work [23], the results of 
this work agree with that of the previous work at least qualitatively. One can 
note that the dynamics of such a system can be realized experimentally. One 




Figure 1: Schematic diagram for initially coiled polymer chain in a locally 
charged channel. The fluidic channel is subjected to an external periodic red 
force A (t) (AC field). 



makes negatively charged coiled polymer chain, then put the polymer within 
positively and negatively charged fluidic channel as shown in Fig. 1. The 
fluidic channel is subjected to an external periodic force (AC field). Since the 
polymer is negatively charged, it encounters a difficulty of crossing through 



the negatively charged part of the channel. Assisted by the thermal back- 
ground kicks along with its conformational change, the polymer ultimately 
overcomes the barrier. The presence of time varying force, may further en- 
hance the rate of crossing. One can sort or manipulate polymer of specific 
coupling constant or chain length by tuning the angular frequency. 

The rest of the paper is organized as follows: in Section 2, we present 
the model. In Section 3, we study the dependence of the crossing rate on 
the model parameters. In section 4, we discuss how the SNR for the flexible 
chain behaves as a function of the model parameters. Section 5 deals with 
summary and conclusion. 

2 The Model 

Let us now consider a flexible polymer chain of size iV which undergoes a 
Brownian motion in one dimensional piecewise linear bistable potential as 
shown in Fig. 1. Considering only nearest-neighbor interaction between 
the polymer segments (the bead spring model), the Langevin equation that 
governs the dynamics of the N beads (n=l,2,3 . . . N) in a highly viscous 
medium under the influence of external potential U(x) is given by 

7-7T = -&(2Xn. - Xn-l ~ X n +l) ^~ + £n(t) (1) 

at ox n 

where the k is the spring (elastic) constant of the chain while 7 denotes the 
friction coefficient. £ n ( t ) is assumed to be Gaussian white noise satisfying 

(£»(*)) = 0, (UQUt + r)) = 2DjS{t) (2) 

with D = ksT is the strength of the thermal noise. The external potential 
each bead experiences is considered to be a piecewise linear potential 



U(x) 



(t/o[=2-l], ifx<-L ; 

U [% + 1], U-L <x<0; 

U [=* + l], ifO<x<L ; {6) 

tfo[£-l], if£>L ; 



where Uq and 2L denote the barrier height and the width of the piecewise 
linear bistable potential, respectively. If one considers only the center of mass 
motion, the second term in Eq. (1) vanishes. For a globular polymer chain, 



where the coupling (spring) constant k is infinity, the Langevin equation (Eq. 
(1)) for the center of mass (cm) takes a simple form 



7V 7 - 






dt 






m, 



if x < —L ; 

if -I/ < a; < 0; 
iVgf+£(*), ifO<x<L ; 
_7Vg + e(i), ifx>L . 



(4) 
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Figure 2: Schematic diagram for initially coiled polymer chain in a piecewise 
linear bistable potential. The potential wells and the barrier top are located 
at x = ±Lo and x — 0, respectively. Due to the thermal background kicks, 
the polymer ultimately crosses the barrier assisted by its conformational 
change along the reaction coordinate. 



3 Polymer escape rate 

We consider a polymer which has initially coiled conformation along the 
piecewise linear bistable potential as shown in Fig. 2. Due to the thermal 
background kicks, the polymer presumably crosses the barrier. The cross- 
ing rate of the chain strictly relies on the chain length, coupling constant, 
barrier height and noise strength. For a compact polymer (k — > oo), the 
monomers feel the same force along the reaction coordinate (see Eq. (4)) 



(3 



and due to lack of chain conformational change at the top of the barrier, 
its escape rate is considerably lower than a polymer of finite k value. On 
the other hand, the jumping rate also depends on the choice of coordinate 
system. The rate for the center of mass motion is a decreasing function of N 
while the escape rate for the portion of polymer segment increases with N. 
The thermally activated rate strictly relies on the chain initial conformation 
along the piecewise linear bistable potential for the case where the micro- 
scopic relaxation time of the polymer is significant compared to the crossing 
time. In this case, initially stretched polymer crosses the barrier faster than 
coiled chain. This is because coiled polymer first stretches before crossing 
the barrier. The degree of stretching depends on the relaxation time of the 
polymer which itself relies on the chain length, coupling constant and the 
thermal background kicks. For short polymer chain surmounting a potential 
barrier that is large compared to D, NUb 3> D, the chain crossing time is 
considerably larger than its relaxation time and hence the crossing time may 
be independent of chains initial conformation. The activated barrier crossing 
considered in this work strictly applicable only to the case where the crossing 
time is much larger than any microscopic relaxation times of the polymer. 

Before we discuss how the flexible polymer chain in the double-well po- 
tential behaves, let us first calculate the crossing rate for a globular polymer 
chain. The problem of surmounting a piecewise linear bistable potential for 
a single Brownian particle in high friction limit was addressed in the work 
[26] . Following the same approach, the mean first passage time MFPT for 
the compact polymer crossing over a high potential barrier is given by 

MFPT=1(^)\^ (5) 

D \NUj K ' 

while the crossing rate 

D,NUo\\=>p (6) 

7 \dlJ y J 
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is the inverse of MFPT. The dependence of the crossing rate r^ or equiva- 
lently the MFPT on the potential height Uq, width of the potential 2Lq and 
chain length N can be analyzed by exploiting Eq. (6). When Uq increases, 
the polymer encounters a difficulty of jumping the piecewise linear bistable 
potential and as a result r^ declines. When N increases, r^ decreases as large 
polymer crosses the potential barrier at the expense of higher thermal kicks. 



On the other hand, as L steps up, the MFPT for the polymer to reach to 
the other side of the well increases which implies r^ monotonously decreases. 

For the chain with a finite coupling constant, we analyze the crossing 
rate via numerical simulation. We introduce dimensionless parameters: x = 
x/L , t = •yL 2 /Uo, k = kL 2 /Uo and i = t/r. Hereafter all the quantities 
are rescaled (dimensionless) so the bars will be dropped. The behavior of 
the system is analyzed by integrating the Langevin equation (1) (employing 
Brownian dynamics simulation). In the simulation, coiled polymer chain 
with N monomers is initially situated in one of the potential wells. Then the 
trajectories for the center of mass of the polymer or the portion of polymer 
segment is simulated by considering different time steps At and time length 
t max . In order to ensure the numerical accuracy, up to 5 x 10 5 ensemble 
averages have been obtained. 

Let us now vary the rescaled spring constant k of the chain. Figure 
3 plots the dependence of r^ on the rescaled k for parameter choice N = 
10, iV = 20, N = 25 and D = 0.5. For small k, since the monomers 
become non- interacting (non-cooperative), r^ tends to be smaller. When 
coupling strength between the monomers k further increases, the tendency 
for interconnected monomers to assist each other increases as a result r^ steps 
up. At certain optimal k optl r^ attains an optimum value. Further increasing 
in k results in a lower crossing rate as larger thermal energy is required to 
drive the rigid chain across the reaction coordinate. The same figure depicts 
that k opt is a decreasing function of N. The coupling constant k opt relies on 
rescaled noise strength D in a manner k = 3DL 2 /l 2 . The crossing rate for 
this entropic chain has been analyzed and compared with the chain of a finite 
k. The numerical analysis shows the rate for entropic chain is considerably 
higher; further details will be reported elsewhere. 

Figure 4a depicts the plot of r^ as a function of N for parameter choice 
D = 0.25, D = 0.6 and k = 0.5. The crossing rate for the center of mass 
motion monotonously decreases with N. On the contrary, the rate of the 
polymer segment r s increases and saturates to a constant value as N and D 
step up as shown in Fig. 4b. 

4 Stochastic resonance 

In the presence of time varying signal, the interplay between noise, sinusoidal 
driving force together with chain flexibility and chain length in the bistable 
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Figure 3: Polymers crossing rate r^ as a function of coupling constant k 
for different chain length N and noise strength D = 0.5. The simulation 
results are obtained directly by integrating Eq. (1). The figure exhibits that 
the rate increases as N declines and attains an optimal value at an optimal 
elastic constant k opt . The optimal coupling constant shifts to the right as N 
decreases. 



system may lead the system into stochastic resonance provided the random 
tracks are adjusted in an optimal way to the recurring external force. Next 
we study the dependence of the SR on the model parameters employing two 
state approximations without considering a continuous diffusion dynamics. 

In the presence of a periodic signal A cos (fit) , the Langevin equation 
that governs the dynamics of the system is given by 



7- 



~dt 



-k[ZX n X n ~\ ^n+l; 



dU{x r 

dx n 



A cos(Qt) + £ n (t) (7) 



where A Q and Q are the amplitude and angular frequency, respectively. 
Employing two state model approach [121 125], two discrete states x(t) 
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Figure 4: (a) Dimensionless rate r^ versus N. (b) The escape rate for polymer 
segment r s versus N. The data are obtained numerically for different rescaled 
noise strength D and coupling constant k = 0.5. 



±L are considered. Let us denote n + and n_ to be the probability to find 
the polymer segment in the right (L ) and in the left (— L ) sides of the 
potential wells, respectively. In the presence time varying signal, the master 
equation that governs the time evolution of n± is given by 



h±(t) = -W±(t)n± + W T {t)n^ 



(8) 



where W + (t) and W_(t) corresponds to the time dependent transition prob- 
ability towards the right (L ) and the left (— L ) sides of the potential wells. 
The time dependant rate [HI [22] takes a simple form 



W- 



± 



r^exp 



UN 



U n D 



A cos (Qt) 



(9) 



where r^ is the Kramers rate for the polymer in the absence of periodic 
force Aq = 0. We consider the case where N ^> D. For sufficiently small 
amplitude, one finds the signal to noise ratio 



SNR = nr k ( 



A NL 
U Q D 



CK2 



(10) 
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Figure 5: (Color online)The dependence of SNRx 10 -6 (for compact polymer 
chain) on N and D for parameter choice Aq = 0.05. The SNR is obtained 
via Eqs. (6) and (10). 

Before we explore how the SNR for globular polymer behaves as a func- 
tion of N, we introduce additional dimensionless parameter: A = AqL/Uq. 
From now on for convenience the bar will be dropped. Figure 5 depicts the 
plot for the SNR for a globular polymer chain versus D and N for fixed 
values of the driving frequency and potential width. The SNR exhibits non- 
montonous noise strength dependence revealing a peak at an optimal noise 
strength D opt . D opt steps up as N inclines. This is because longer polymer 
crosses the barrier at the expense of higher thermal excitation. The same 
figure exhibits that the SNR peaks at optimum chain length N opt . This 
suggests novel means of manipulating not only for biopolymers, but also for 
proteins (vesicles) of different size. One can note that since compact polymer 
lacks the flexibility to conform to the piecewise linear bistable potential, the 
peak of the SNR is less than that of a flexible polymer chain. 

For the chain with a finite coupling constant k, the resonance behavior of 
the system is examined numerically. Figure 6 shows the plot for the SNRk 
as a function of k for different values of N. Other parameters are fixed as 
D = 0.25 and A = 0.1. For small k, the monomers tend to be non-interacting 
as a result the SNR becomes very small. The SNR peaks at an optimal chain 
length kopt and further increasing in k leads again to a smaller SNR as the 
rigid polymer is not flexible enough to adjust itself with the time varying 



11 



SNR, 



0.00006 




^e*H 


0.00004 


4 




0.00003 


/ # 

T f 

1 


\\ 


0.00001 



j i 

f 


X 



0.2 0.4 0.6 0.8 1 1.2 1.4 

k 

Figure 6: The dependence of SNRk as a function of coupling constant k 
for different chain length N, for noise strength D = 0.25 and amplitude 
Aq = 0.1. The simulation results are obtained directly by integrating Eq. 
(1). The figure exhibits that the SNR increases as N declines and attains an 
optimal value at an optimal elastic constant k opt . Unlike the corresponding 
rate, the optimal coupling constant shifts to the left as N decreases. 



force. The optimal spring constant k opt increases as N increases. 

Figure 7a plots the rate SNRk versus N. We observe an increase in the 
signal to noise ratio when the number of monomers increases. The SNR 
peaks at an optimal chain length N opt and declines again as the chain length 
increases. On the contrary, SNR S monotonously increases with TV (see Fig. 
7b). It is worth noting that since the flexible polymer chain responds in 
more cooperative and coherent manner to the time varying force, it exhibits 
a higher resonance peak than a globular polymer. 

Finally we emphasize that even though the model and its numerical ap- 
proach is completely different from our previous work [23], the results agree 
with that of the previous work at least qualitatively. For instance, similar 
to this work, in our previous work we found that the response of signal to 
the background noise strength is significant at optimum values of D opt , K opt 
and Nopt. Our previous work [23] dealt with the study of SR for a flexible 
polymer chain on Kramers potential. Depending on the chain length iV and 
spring constant k, the chain takes either coiled or stretched conformation at 
the top of the barrier. When the chain is either in coiled or stretched state, 
the resonance is much larger than globular state revealing the intrinsic flex- 
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Figure 7: (a) SNRk versus N. (b) SNR S for polymer segment as a function of 
N. The data are obtained numerically for noise strength D = 0.5, amplitude 
A = 0.1 and coupling constant k = 0.5. 

ibility of the chain facilitating faster crossing. In the globular limit k — > oo, 
the resonance becomes small as the compact polymer lacks the flexibility to 
conform to the driving force. Furthermore, we showed that the power am- 
plification peaks at an optimal chain length and elastic constant as well as 
at an optimal noise strength. 

5 Summary and conclusion 

In summary, in this work we explore the crossing rate and stochastic res- 
onance of a flexible polymer chain in a piecewise linear bistable potential. 
We investigate the chain's escape rate as a function of different model pa- 
rameters. Tfc monotonously decreases with N while the rate r s increases and 
saturates to a constant rate when N steps up. Due to lack of cooperation 
between the monomers, the rate is considerably small for a smaller k. The 
crossing rate increases as k goes up. At certain k opt , r^ attains an optimal 
value. Further increasing in k results in a lower rate as rigid polymer crosses 
the barrier at the expense of higher thermal kicks. 

Employing two state approximation, the dependence of SNR as a func- 
tion of model parameters is studied. For a globular polymer chain we show 
the response of signal to the background noise strength is significant at op- 
timum values of D opt and N opt which suggests novel means of manipulating 
proteins (vesicles) of different size. On the other hand, for the chain with a 
finite coupling constant k, the SNR exhibits an optimal value at an optimal 
kopt- The signal to noise ratio for the center of mass SNRk exhibits an op- 
timal value at an optimum N opt . The SNR for the polymer segment SNR S 
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monotonously increases with N. 

In conclusion, in this work, by introducing a different model system than 
the previous work [23], we recapture the previous results at least qualita- 
tively. Not only we assess the resonance property of the system, but we also 
further explore the barrier crossing rate by varying different model param- 
eters. Since polymers are interconnected and flexible systems, they exhibit 
interesting cooperative dynamics when they are exposed to time varying ex- 
ternal fields and noises. Understanding their dynamics is crucial not only for 
novel means of manipulating proteins such as DNA or RNA molecules in a 
nanofludic or microfiudic channels, but also to understand how such systems 
self-organize their flexible degrees of freedom. Thus, this theoretical study is 
crucial not only for the fundamental understanding of polymer physics, but 
also provides a basic paradigm in which to understand the self-organization 
and cooperativity induced by the chain flexibility and fluctuations. 
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